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In the category of motions preserving the angular momentum's direction, Gorringe and Leach 
exhibited two classes of differential equations having elliptical orbits. After enlarging slightly these 
classes, we show that they are related by a duality correspondence of the Arnold- Vassiliev type. The 
specific associated conserved quantities (Laplace-Runge-Lenz vector and Fradkin-Jauch-Hill tensor) 
are then dual reflections one of the other. 



PACS numbers: 

I. INTRODUCTION 

In 1993, Gorringe and Leach [ij, exhibited two classes of differential equations incorporating drag terms which 
have closed elliptical orbits, generalizing then previous results of Jezewski and Mittleman and Leach [1]. Both 

possess conserved quantities which extend the Laplace-Runge-Lenz vector and Fradkin-Jauch-Hill tensor respectively. 
These two classes belong to a broader category ofplanar motions incorporating velocity dependent terms submitted 
to certain constraint on their coefRcients [1, @, 0, Q ■ 

In this paper we show that the above classes can be slightly enlarged. To these generalized Gorringe-Leach equations 
are associated two conserved quantities : a pseudo-energy and a pseudo- angular momentum. We obtain two types 
(H and K) of generalized Gorringe-Leach equations presenting closed orbits for every values of these quantities. In 
the spherically symmetrical case, we obtain compact analytical formulas for the periods which, if restricted to the 
standard case, recover the results of ref. The H and K types belongs to a larger category of generalized Gorringe- 
Leach equations possessing duality properties (in the Arnold- Vassiliev sense [ol. Iioj). In this category, the equations 
can be gathered in classes indexed by a characteristic real parameter ly. Each i/-class possesses an associated dual 
/i-class, with /i — —j^. H type and K type generalized Gorringe-Leach equations are then shown to be dual of each 

other. As in the conservative case [l^l, the pseudo Laplace-Runge-Lenz vector associated to the K type equations is 
also proportional to the dual transform of the pseudo Fradkin-Jauch-Hill tensor of the H type equations. 

II. COMPLEX DESCRIPTION OF MOTIONS WITH CONSERVATION OF THE DIRECTION OF 

ANGULAR MOMENTUM 

As was shown in [s!], the most general form for the equation of a 3-dimensional motion ~r{t) for which the angular 
momentum L = L~e l conserves its direction is given by 

~r + h~r + gT' = 0, (1) 

where h and g are two arbitrary scalars depending on time t. 

As the motion is confined to the plane orthogonal to ~f?L, we can adopt a complex representation for the position 

The functions h and g are then represented by two real arbitrary functions of t through z{t) and z{t), not necessarily 
analytical (that is potentially dependent on z(t) and z(<)). 

We apply our interest more specifically to the autonomous case {g and h do not depend explicitly on t), to which 
we add the constraint that the last term is a function of z and z only. We then arrive to the following equation for 
our planar motion z[t) : 

z + h{z,z, z^'z)z + g{z,l)z — (2) 



g and h being two arbitrary real-valued functions. 
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III. EULER-SUNDMAN REPARAMETRIZATION 

We now perform an Eulcr-Sundman reparametrization t G — > s £ M+ (where the correspondence is one to one 
and increasing) of our motion, z(t). We put 

s = s{t,z{t),z{t)) (3) 

with 

If we write ^ — f\ Eq. [5] becomes 

,, s + sh(z,z, sz' , sz') , ( 
z" H ■ -z' + - 

(^) 



z"+ ^ + ^^("'y'^"\ -+gi^z = (5) 



s 



in which s is called the pseudotime and z{s) the pseudomotion. 

If we choose the reparametrization in such a way that 

s + sh{z,z, sz' , s'z') — 0, (6) 

then Eq. [5] for the pseudomotion becomes 

z" + g{z,z) (e^ ^ h(z,z,z^z)dt\ 

V / t=t(s) 

The problem is considerably simplified if we restrict /i to be a total derivative : 

h{z,z, z,^) = H{z,z) (8) 
In this particular case the initial Eql2]is written as 

z + H{z,z)z + g{z,z)z = (9) 

The reparametrization is given by : 

ds =e-"'-'-^Ut (10) 
and the pseudomotion Eq. [7l takes the following simple form 

z" + 5(z,z)e2^("^^)z = (11) 

which is the equation of an autonomous conservative motion. The term g{z,z)e'^^^^''^'^z is called the pseudoforce. 
From now we always place ourselves in this case. 

IV. RADIAL PSEUDOFORCE 

Consider the case in which the pseudoforce, g{z,z)e^-^^^''^^z, is derived from a radial real-valued potential U{r): 

VC/(r) = = (p{r)z = g{z,z)e^"^''~^z (12) 

where 

/ N 1 9U(r) 
r or 

The initial motion Eq. [9l takes the form: 

'z + H{z,z)z + ip{r)e-^"^^'^'>z = (14) 
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which we call the generalized Gorringe-Leach equation fl*]. 

To the pseudomotion wc can associate a pseudo-angular momentum, 



and a pseudoenergy, 



£ = X 7^' = (z'z — zz') ~e (15) 
2i 



£ = i|yf + [/(r)^i|z'|Vt/(r)), (16) 



which are constants of the pseudomotion. 

Consquently the original motion possesses also the two conserved quantities: 



1 £: = e2^(-'^)^ + C/(r). 

V. BERTRAND'S THEOREM 

Applying Bertrand's theorem [HI, [13] to the pseudomotion, we deduce immediately that the only pseudomotions of 
which orbits are closed for every value of the characteristic parameters (given by the pseudo-energy and the pseudo- 
angular momentum), are those associated to the Hooke and Kepler potentials, 

UH{r) = hr-" (18) 

and 

UK(r) = -- (19) 
r 

In both cases the orbits are ellipses, centered in O in the H case and having a focus in O in the K case. 

The orbits of the pseudomotion being also those of the initial motion (the two motions differing only by a 
reparametrization), we deduce that the only generalized Gorringe-Leach Eq |141 which always (that is for every value 
of the conserved pseudo-energy £ and pseudo-angular momentum C) admit closed orbits, are those for which g{z,'z) 
is of the form: 



gH 



(z,z) = fce-2^("'") (20) 



or 



gK(z,z) = ^e-^^(--). (21) 
In other words, the initial equation must be of the form 

z + H{z^)z + ke^'^^'-^-^h = 0, type H equation, (22) 

or: 

z + ff(^)z + fcr~^e"^^(''^^)z = 0, type K equation. (23) 

VI. SPHERICALLY SYMMETRICAL CASE 

In the spherically symmetrical case H and g depend only on the radial variable, r = -s/z^, ie 



H{z,z) = H{zz) = a{r) 
9{z,z) = g{zz) = u{r), 



(24) 
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and we are in the case in which the pseudoforce derives from a radial potential given by 

U{r) = J ruir)e^°'^''Ur. (25) 

In this case the initial EqUHis written as 

z + a{r)z + u{r)z = 0. (26) 

The two conserved quantities Eqs. ?? become: 

I IP (27) 

and the type H, and K Eqs. I22l23l takc the forms 

z + a(r)z + fce^^"^^-'z = 0, type H equation, , , 

z + a{r)z + kr^^e^^"''''^^ z — 0, type K equation. 

Note that, if we restrict ourselves to the case of a logarithmic function for a(r), we recover here in a very direct 
manner the original Gorringe-Leach results P, . 

VII. ORBITAL PERIOD IN THE SPHERICAL SYMMETRICAL CASE 

The radial component Ws) of the pseudomotion is given by integration of Barrow's differential formula applied to 
the radial pseudomotion |12l | , which gives 

ds 1 e"('')dr 
dt^ — = ^ , 29 



where the radial effective potential Vl' (r) is 



Vcir) = U{r) + ^. (30) 
2r^ 



When the orbit is bounded, the radial variable r oscillates between the extremal values a et & (pericentral and 
apocentral radii) which are roots of the numerical equation: 

£-Vc{r) = Q. (31) 

The orbital period is 

We are always in this case for the type H and K equations. 

If, for the corresponding elliptical orbits, we note A the minor axis and B the major axis, we have : 



* In the H case {£ > 0) V^{r) = \kr^ + 



A = a 
B = b 



(33) 



and 



2 f'' re"('')dr , , 

Th = ^ , 34 
VkJa ^(r2 - a2) (62 _ r^) 



* In the K case {£ < 0) Vf{r) = -f ' ^ 



b = B + y/ W^^ 

a = B- VW~A^ 



(35) 
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and 



T^K = -1= / =■ 36 



By a straightforward changes of variables we finally obtain : 

Th = ^ x-i (1 - xy^ exp (a(V(^'2-a2)a; + a2)) da; 
= ^^5^ -^o' (1 ~ 2;)"^ (a; + _a_) (a((6 - a) a; + a)) dx 

VIII. GORRINGE-LEACH EQUATIONS 

Gorringe-Leach equations [l| are equations of the type considered with 



The initial equations of motion presenting closed orbits are then (see Eq. ??) 



and 



As for the corresponding orbital periods Eq. [37l they are written in these cases as 

ra^t / a 1 



and 



(37) 



a(r) = -|ln(r). (38) 



z- ^-z + A:r"z = 0, H type, (39) 



'z z + kr^-^z = 0, K type. (40) 

2 r 



Th = ^^f(-,-,1, ^ (41) 



„ \/2Tra^ 2 /a 1 a — b\ 

Tk = -^-^^^F[ , (42) 



\2 '2' a 

where F is an hypergeometric function. 

Note that the hypergeometric functions of the form F {v, f3, 2/3, z) can be rewritten in terms of Legendre functions 
via : 

Fi,.,P,2p,z) = 2"'-'Ti(3+^)z^-^{l ~ z)^Fj_:/_, (t^) ' ^^^^ 
Under this form we recover the results of Gorringe and Leach p| , namely 

Th = ^(1 - z)-fp., ( i=l') (44) 



and 
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IX. ARNOLD- VASSILIEV DUALITY 

In order for the pseudomotion potential be an Arnold- Vassiliev potential and then be dualizable in the 

Arnold- Vassiliev sense, it has to be of the type 

U{z,z) ^ k\u{z)\^ eR (46) 
It is easy to show that, if we suppose u{z) analytical, this necessitates 

z^^ = cste^UR, (47) 

that is, U must be a central power law potential namely 

U{z,z) = Ar'' (48) 

and the pseudomotion equation is (k = vA) 

z" + kr^-^z = 0. (49) 

In other words the generalized Gorringe-Leach equations (Eq. [T4|l which are dualizable in the Arnold- Vassiliev 
sense have the form 

'z -f H{z,z)z + /jr'^-2e-2-f^(^'")z = 0. (50) 

Such an equation is said to be of class v. In this nomenclature an H-type equation is then an equation of class 2 and 
a K-type equation an equation of class -1. 

The pseudomotion corresponding to an equation of class v admits a dual in the Arnold- Vassiliev sense, the associated 
potential of which is given by 

V{w,w) ^ Bpf", (51) 

where p^\w\, B = -f / (l f )' , (l + f ) (l + f ) - 1. 
The dual pseudomotion equation is then 



w" + Kr^-^w = (52) 

with K = nB = -pE' / (1 + f )^ 

The correspondence between the position variables is 

^ (l + I) zi+f (53) 

and the dual pseudotimes are related by the Euler-Sundman reparametrization 

cls={l + ^ ^ p'^da. (54) 

To the pseudo- motion, w{a), is also associated a "true" motion, w(r), satisfying a generalized Gorringe-Leach 
equation of the form 

w + H{z,z)w + Kr^~2g-2H(2,2)^ ^ (-55^ 

H{z,z) being an arbitrary real valued function. 

The dot indicates here the derivative with respect to t, where the "true" time, r, is related to the pseudotime, cr, 
by the Euler-Sundman reparametrization 



da ^e-"^'-^'-^dT. 



(56) 
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This equation is an equation of class fj,. We can then see that the generahzed Gorringe-Leach equations, if not strictly 
dualizable, can, however, be grouped into classes indexed by a real characteristic exponent, the v class being in dual 
correspondence with the class ^, such that : 



2/ V 2 

The classes 2 and —1 are dual in the sense that they are linked by a Levi-Civita change of coordinates, this 
preserving the elliptical structure of the orbits. The fact that the equations of H and K types possess the same orbital 
characteristics is not fortuitous, but is a direct consequence of this duality. 

As was demonstrated in Td] , in complex formulation the existence of an additional conserved quantity for a motion 
in the Hooke potential (the Fradkin-Jauch-Hill tensor) is evident. If we apply this result to a type H pseudomotion, 
we deduce immediately that a H type generalized Gorringe-Leach equation admits besides its pseudo-energy, £, an 
additional conserved complex quantity, namely 

r = 1 (z'f + \kz^ = e2^(--)i (z)^ + \kz\ (58) 
(The components of the 2D associated pseudo-FJH tensor are given by T^ja; — Re(£ -|-T), Tyy — Re (£ — T) and 

T,j,y = Im(T) =Ty^.) 

Following [To| the dual K-type generalized Gorringe-Leach equation iii + H{w,w)'w + kp^^e^^^^^''"''w = 0, admits 
in addition to its pseudo-energy £, the following conserved complex quantity : 

A=L^ - !^ = e^^(--) ^^w'L - ^ (59) 

2k k P k P 

which is the image of T via the dual transformation (EqslSSland lM)) . with 8 — — ^ and k = We clearly recognize 
in A the complex formulation of the pseudo-Laplace-Runge-Lenz vector pj] : 

A = ^L xy (60) 
k P 
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